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Feedback stabilization for the mass balance 
equations of a food extrusion process 

Mamadou DIAGNE Peipei SHANG Zhiqiang WANG 


Abstract —In this paper, we study the stabilization problem 
for a food extrusion process in the isothermal case. The model 
expresses the mass conservation in the extruder chamber and 
consists of a hyperbolic Partial Differential Equation (PDE) 
and a nonlinear Ordinary Differential Equation (ODE) whose 
dynamics describes the evolution of a moving interface. By using 
a Lyapunov approach, we obtain the exponential stabilization for 
the closed-loop system under natural feedback controls through 
indirect measurements. 

Index Terms —Feedback stabilization, hyperbolic system, mov¬ 
ing interface, Lyapunov approach. 

I. Introduction 

CREW extruders have become very popular for their abil¬ 
ity to manufacture food and plastics products with desired 
shapes and properties. Due to the strong interaction between 
the mass, the energy and the momentum balances occurring in 
those processes, the design of efficient controllers still remains 
a hard task at the industrial level. So far, the control oriented 
model of extruders are issued from some black box model 
of limited operational validity. Following the objectives of 
the control, these models describes the extruder’s temperature 
and flow rate at the die output or the pressure dynamics 
based on single input and single output or multiple input 
and multiple output system. Generally, extrusion processes are 
controlled using PID (TO) , p~3) , | ~j~8| or predictive controllers 
(20) with oversimplified or empirical models. In (13) , the 
volumetric expansion of the extrudate correlated to the die 
temperature and pressure and the specific mechanical energy 
is chosen as the key product quality to be controlled. The 
authors study the performance of the PI controller based on the 
regulation of the die pressure using feed rate as a manipulative 
variable and show that the response of an improperly tuned 
controller may be too sluggish on one hand, or too oscillatory 
on the other hand. First-order, second-order and Lead-lag 
Laplace transfer-function are exploited in HD to design a 
feedforward controller for a twin-screw food extrusion process 
to reduce the effect of feed rate and feed moisture con¬ 
tent variations on the die pressure. The experimental results 
showed that the die pressure response varies with the operation 
conditions, so a single gain value would not be suitable 
for all operating conditions. The complexity of the process 
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suggests adaptive control methods because of the frequent 
changes occurring during the extruder operation in the feed 
composition E3- ED uses second-order transfer functions 
while emphasizing the difficulty in implementing those types 
of model-based controllers due to the strong influence of 
all the manipulated inputs and measurable process variables. 
Therefore, very intelligent controllers need to be constructed 
for extrusion cooking process based on a control algorithm 
developed from process experience. We mention that the 
transport delays, the strong interactions and the non-linearities 
make it difficult to control such systems with PID controllers. 
Predictive controllers might offer better performances but are 
somewhat difficult to implement ea¬ 
rn the present work, we consider the stabilization of the die 
pressure to the desired setpoint in a food extrusion process. 
The controller is constructed based on a bi-zone model derived 
from the computation of the conservation of mass in the 
extruder under the assumption of constant temperature and 
viscosity. A Geometric decomposition of the extruder length 
in Partially Filled Zone (PFZ) and Fully Filled Zone (FFZ) 
allows to describe the process by a transport equation and 
a pressure-gradient equation defined on complementary time 
varying spatial domains. The domains are coupled by a moving 
interface whose dynamics is governed by an ODE representing 
the total mass balance in an extruder. We emphasize that 
an accurate die pressure regulation is critical to ensure the 
uniformity of the extruded melt and is strongly related to the 
quality of the final product. We propose a suitable feedback 
control laws together with practical measurements as output 
such that the solution of the closed-loop system converges to 
a desired steady-state or equilibrium asymptotically. 

The stabilization problems for hyperbolic systems has been 
widely studied in the literature. The first approach relies on 
careful analysis of the classical solutions along the charac¬ 
teristics. We refer to Greenberg and Li ED in the case of 
second-order system of conservation laws and more general 
situations on nth-order systems by Li in (16) . 

Another approach based on Lyapunov techniques was intro¬ 
duced by Coron et al. in [4|. This approach was improved by 
in (3) where a strict Lyapunov function in terms of Riemann 
invariants was constructed and its time derivative can be 
made negative definite by choosing properly the boundary 
conditions. The Lyapunov function is very useful to analyze 
nonlinear hyperbolic systems of conservation laws because 
of its rubustness, see Q-@, ©. ©, (22), (23) for a wide 
range of applications to various models. Among which, we are 
interested in a physical model for the extrusion process which 
occurs very often in polymer material and food production. 
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The main contribution of this paper is to establish the 
exponential stabilization for the extrusion model under natural 
feedbacks by a Lyapunov function approach motivated by [3j, 
(23). The difficulties come from three aspects: 1) The 
domains on which the conservation laws are defined depend on 
the solution through the dynamical interface; 2) The nonlinear 
coupling of the dynamics of the interface and the filling ratio 
in the Partially Filled Zone (PFZ) is not standard; 3) The 
measurements and the feedbacks are natural, however, some 
of the measurements are indirect, i.e., the measurements are 
not a part of solution but given through a indirect relation of 
solution. 

The organization of the paper is as follows. In Section [II] 
we give the physical description of the model. The main result 
on stabilization (Theorem [TJ and its proof are given in Section 
[III] Numerical simulations are provided in Section [IV] Finally 
in Section [V] we give our conclusion and some perspectives 
in control of the extrusion process. 

II. Description of the model 


convection speed in the two regions. The flow rate in the 
FFZ is constant and equal to the die flow rate F d {t) which is 
proportional to the pressure difference A P{t) := P(t. L) — P 0 
where Pq denotes the atmospheric pressure. For more detailed 
physical description of the model and definition of all the 
parameters appeared below, one can refer to 0> §• 

In this work, the stabilization of ( l[t ), f p (t , x)) with the help 
of the actuated screw speed N(t) and inlet flow rate Fi n (t) 
is established based on feedbacks that depend on the pressure 
difference A P(t) that is a practically useful measurement for 
the system. Considering the following change of variables 


x>y = 


x 

W) 


in PFZ and x >-> y 


x — l{t) 
L — l(t) 


in FFZ , (1) 


respectively, the time varying domains ([0, l(t)], [l(t),L]) can 
be transformed to the fixed domain [0,1] in space. For the sake 
of simplicity, we still denote by x the space variable instead of 
y. More precisely, we consider the stabilization problem for 
the corresponding normalized system on the spatial domain 
[0,1], The interface dynamics which arises from a total mass 
balance writes 


An extruder is a process used for manufacturing objects 
with fixed shapes and specific properties, see Fig. [I] One 
or two Archimedean screws are rotating inside the barrel in 
order to convect the extruded material from the feed to the 
die exit. In food or polymers extrusion processes, the ultimate 
control systems involved manipulation of screw speed, feed 
rate, inlet solvent fraction, and barrel temperature for the 
regulation of moisture content, temperature and viscosity of 
the finite product, residence time and die flow, etc. In this 


feed moving heater 



Fig. 1. A single-screw screw extruder 


paper, we consider the mass balances model 0,(8) motivated 
by (14), (T5) for cooking extrusion process. In this case, the 
material convection along the extruder chamber of length L 
is described in two zones: the PFZ ([0,((t)] in space) and 
a FFZ ([l(t),L\ in space) separated by a moving interface 
l(t). The PFZ and the FFZ appear due to the die resistance 
that provokes an accumulation phenomena and high pressure 
need to be built-up to evict the extrudate out of the die. By 
the mass conservation principle the convection in the PFZ is 
described by the evolution of the filling ratio f p (t,x) for an 
homogeneous melt. The melt convection speed in the PFZ , 
namely, a p depends on the screw speed N(t) whereas the 
FFZ transport velocity is related to the die pressure P(t,L ): 
a p / Of. Under the assumption of constant viscosity ?/ 
along the extruder, the dynamics of the moving interface l(t) 
is governed by an ODE arising from the difference of the 


lit) = HKt), N(t),f p (t, 1)), in K+ = (0, oo), 

m = i °, 


where 


F(l(t),N(t)J p (t,l)) = 


^AP(f) - p 0 V eff N(t)f p (t, 1) 


Po^effi 1 fp{t,l)) 


(3) 


Kd , pa, Tj denote the die conductance, the melt density and the 
viscosity, respectively. V e ff and S e f f are the effective volume 
and section between a screw element and the extruder barrel, 
respectively. Assuming a constant viscosity along the extruder 
(the isothermal case), the relation 


A P{t)=P{N{t)At)):= 


r)PoV e ffN(t)(L - l(t)) 
Bp 0 + K d (L — l(t)) 


(4) 


is determined by integrating the pressure-gradient equation 
corresponding to the momentum balance in the FFZ and 
considering a pressure continuity coupling relation at the 
normalized interface, namely, P(0,i) = Pa in the PFZ [81. 
The filling ratio in the PFZ writes 


f d t f p (t,x) + a p d x f p {t,x ) = 0, in K+ x (0,1), 
lf P (0,x) = f°(x), in (0,1), 

(f P (t,0) = Po vT t %(tV inR+ ’ 

where 

_ c N(t)-xPm,N(t)J p (t,l)) 


(5) 


( 6 ) 


III. Main result and its proof 

Let us define the constant equilibrium l e £ (0 ,L), / Pe £ 
(0,1), N e , A P e , a Ps and F ine by F(l e ,N e ,f p J = 0, 
A P e = V(l e ,Ne), a P e = ^,Fin e = p 0 V e ffN e f p _ s . Denote 
the difference ((f) := l(t) — l e , N{t) := N(t) — N e ,f p (t,x) := 
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f P (t,x) - f Pc ,F in (t ) := F in (t) - F inB ,AP(t) := A P(t) - 

A P e and the constants 


f(ai , 02 , 03 )= 


(le,N e ,f Pe )’ 


(7) 


The linear feedback law that we use is the following one: 


N(t) = k, • A P(t), F in (t ) = fc 2 • A P(t), ( 8 ) 


where A P(t), thus A P(t), is measurable. The aim of stabi¬ 
lization is to find constants k \, k 2 £ R such that the closed- 
loop system 0 and 0 with feedback 0 is asymptotically 
stable, i.e., ( l(t),f p (t , •)) —► 0 as t —» 00 . 

Concerning the well-posedness of the Cauchy problem ([2]) 
and <[5j with feedback ([ 8 }, we have the following proposition. 


Proposition 1. Let k \, k 2 £ R be fixed. There exists e > 0 
such that for any l 0 £ R, / p ° £ H 2 ( 0,1) satisfying \l° — l e \ 2 + 
ll/p(') — fpJ\ 2 H 2 (o l) ^ £ ’ an d ^ le compatibility conditions 
at the point ( t,x ) = ( 0 , 0 ), system 0 and 0 with 0 has a 
unique solution (l,f p ) £ W 1 ,oo ([0, T)) x C°([0, T); H 2 (0, 1)) 
for some T £ (0,oo]. Moreover, if \l(t) — l e \ 2 + \\f P (t, •) — 
/pj|^ 2 ( 0>1 ) < efor all t £ [0,T), then T = 00 . 


Remark 1. The compatibility conditions at the point (t, x) = 
( 0 , 0 ) are the following: 


= Fin ( 0 ) 

Jp[) PoV eff N{0y 
F in (0)N(0) ~ F in (0)N(0) 
Po V eff N{ 0) 


CN{p) 

l 0 


(9) 


fp ( 0 ) = 0 , ( 10 ) 


where N(0), Fi n (0) are determined by (JT]», | 8 ]( with 1(0) = 1°, 
while N( 0), Fi n (0) are determined by differentiating 0 and 
<[0» together with 1(0) = F(l°, N(0), fp (1)). 


The proof of Proposition [T] is based on fixed point argument 
and one can refer to [ 8 j for the well-posedness of the corre¬ 
sponding open-loop system. Our main result on stabilization 
of the interface position l(t) and the filling ratio f p (t,x) is 
the following theorem. 


Theorem 1. Suppose that there exist k \, k 2 £ R such that 


CLl + 


fcia 2 &i 
1 - kib 2 


< 0 , 


«3^l(fc2 ~ fp c PoVeffkl) 
PoV e ffN e ( 1 - kib 2 ) 


< 


(L\ + 


kici 2 bi 
1 - kib 2 


( 11 ) 

( 12 ) 


where ai, a 2 , 03 , bi, b 2 are given in 0 . Then, the nonlinear 
system 0 and 0 is locally exponentially stable under the 
feedback 0 , i.e., there exist constants e > 0, M > 0 and 
to > 0 such that for any 1° £ R, /° £ H 2 (0 ,1) satisfying 


\l°~ (e| 2 + !!/“(•) -/ Pe ll|p(0,l) O, (13) 


and the compatibility conditions at the point (t , x) = ( 0 , 0 ), 
the solution of 0 and 0 with 0 satisfies 

\l(t) — l e |“ + || fp(t, •) - /pj|tf2( 0 ,l) 

< Me-^{\1° -l e | 2 + ||/°(-) - / p J| 2 ff 2 (0il) ), Vf ^ 0. 

(14) 


Before the proof of Theorem [I] we give several remarks. 

Remark 2 . The existence of k \, k 2 £ R satisfying >[TT])-([T 2 |) 
depends on ai, a 2 , 03 , b±, b 2 , see Proposition^below. Its proof 
is given in Appendix. 

Proposition 2. There exist k\. k 2 £ R such that (| 11 1-(| 1 2[i hold 
if and only if either 1) a± < 0 ; or 2) a\ ^ 0 and a 2 b\ ^ 0 . 

Remark 3. The measurement on A P(t) is of practical reason, 
thus the feedback 0 is indirect in the sense that the measure¬ 
ments are made not on the solution (l(t), f p (t, x)) itself. 

Remark 4. The proof of Theorem [7] relies on a Lyapunov func¬ 
tion approach. The weight as e~ liX is essential to get a strict 
Lyapunov function. One can refer to the stabilization results 
by such weighted Lyapunov functions, for quite general linear 
hyperbolic systems in for one dimensional Euler 

equation with variable coefficients in (00 ; for a conservation 
law with nonlocal velocity in /HP- 

Proof of Theorem [lj The construction of the Lyapunov 
functions is divided into three steps. 

Step 1. The stabilization of l(t) and f p (t,-) in L 2 (0,1). 

Let 

V 0 (t)=P(t), V 1 (t)= f e-^ x f 2 p (t,x)dx, (15) 

J 0 

where ji > 0 is a constant to be chosen later. 

Lemma 1. There exist positive constants A\, 71 , /3q, /3i, <5i 
such that the following estimates hold for every solution to 
system 0 and 0 with 0 

Vo(i) + A\V\(t) ^ — (A) + o(l))Vo(t) — (At + o(l))Ci(t) 

— (A + o(l))fp(t, 1 ), (16) 

where o(l) represents various terms which tend to 0 when 

\m,N(t)j p (t,i))\^o. 

Proof of Lemma |T} By definition of the equilibrium 

(l e ,N e , f p f) and the constants (ai, a 2 , ( 13 , b±, b 2 ), it is easy 
to get by expansion that 


F(l(t),N(t), f p (t , 1)) = (ai + o(l)) l(t) + (a 2 + o(l))N(t) 



+ (a 3 + 


(17) 

AP(f) — (&i + o(l)) l(t) + (b 2 + 

o(l)) N(t). 

(18) 

Furthermore, it follows 

from <[ 8 ]> and (fl 8 | that 


A P = 

( h 

Vl — kib 2 

+ o(l))/(f), 

(19) 

N(t) = 

( hh 
vl — k\b 2 

+ °(l))U/ : )) 

( 20 ) 

Fin (t) = 

( k 2 bi 

VI — k\b 2 

+ Al))^)- 

( 21 ) 


Differentiating Vo (t) with respect to t and using Q, ( fT7| ) 
and ( [ 20 ] >, one easily gets that 

V 0 (t) =2l(t)-F(l(t),N(t),f p (t,l)) 

=(20 o + o(l))P(t) + (2a 3 + o(l))I(t)f p (t, 1). (22) 

, . k!a 2 bi 

where 0 O = at H- r~r~ < 0 - 

1 - kfij 2 
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On the other hand, © and {17} yield that 


QZp CT Pe T o(l), 


= 0 ( 1 ). 


V 1 (t)=BT 1 + [ (-71 a p + a p Je llX f p (t,x) 
Jo 

=BTi-(7iap e +o(l))Vi(t), 


da; 


where 


BTl = —e~ 71 (a Pe + o(l))j?(i, 1) + (a Pe + o(l))j?(i, 0). 


Note that by {5},{20}-(|2Tj), we have 


f P (t, 0 ) = 


F in (t) 


PoV e ffN(t) 


-f P ' = (0i+o{i)m, 


+ (a Pe 6 1 + o(l))Vo(t). 


29o 4 - Aia Pe 9 2 0,3 

<13 — A\e~ 11 a p 


fp*t a pfp** + <*pjp. — 0 ; 

f P A°’ x ) = fp( x )» 

f P At,o) = ~ fpt{t,0) 


where 


Let 


fp t (t, °) = 


&p\x=0 

F m (t)N(t) - F in (t)N(t) 
p 0 V eff N 2 (t) 


V 2 (t) < - (ft + o(l))V 2 (i) - (fc + o(l))/£(i, 1) 


+ 9 2 (V 0 (t) + f (t, 1)). 


Then it follows from ([2]), {7}, {T7} , {19} and ( [33] ) that 


(23) 


Differentiating V\ (t) gives, from {TT}, ((23), ((23) and integra¬ 
tion by parts, that 


AP(t) = 


dV 1 

91 1 mMt)) 

1 -k,FP\ 

1 ft-1 atv l(i(t),jv(t)) 


■i(t) 


= 0 (l)(Z(f) + lV(t) + / p (f,l)), 


(34) 


(24) 


(25) 


(26) 


where 0 ( 1 ) denotes various terms which are uniformly 
bounded when \(l(t), N(t), f p (t,l))\ —> 0. 

Combining {70} , l|29|-(|3()[i and l|33[i-(|34ji, we get easily that 

f Px (t,0) (35) 

Differentiating ( (31) results in, by ( [23] ), {23} and ( (29) , that 

V 2 (t) =BT 2 + f (— 72 ^ - a Px )e~ l2X f p ^ (t, x) dx 
Jo 

=BT 2 + (-72 a Ps + o(l))V 2 (t), (36) 

where 


where 01 = ■ Combining ((23), ((24). {25), 

we get consequently 

Vi(t) = - (7i« Pe + o(l))Vi(t) - (e~ 7l a Pe + o(l))/ 2 (f, 1) 


BT 2 = (-£ 


o(l))/^(t,l) 


+ (a Pe + o(l))f Px (t, 0 ). 

Thanks to {23} and ( (35) , ( [37] ) can be rewritten 


as 


(27) 


5T 2 =(—e 


—72 


G(.n 




Under the assumption of (JTT)-([T 2 ), it is easy to get the 
existence of A\ > 0 and 71 > 0 (suitably small) such that 


+ 0 (l)(Z 2 (i) + / 2 (M)) 


(28) 


which ends the proof of Lemma [2] with 

Step 3. The stabilization of f Pa . x (t, ■) in L 2 { 0,1). 


(37) 


(38) 

□ 


is negative definite. This concludes the proof of Lemma [T] with 
{22) and {27}. □ 

Step 2. The stabilization of f Pn .(t , •) in L 2 { 0,1). 
Differentiating © with respect to x, we get 


By differentiating 

fp**t + a pfp*. 


note that a Pxas = 0 , we derive that 


4~ 2 otp x fp xx — 0, 


fp**(°> x ) = fp 0*0 > (39) 

-JWt (C 0) + u P * l*=o f P * (t, 0) 


f P **(t, 0) = 


where 


(29) f Pxt (t, 0) = -f 


d / -1 


|cc =0 


ZP; 

Let 


dt Vet 


p|x =0 


PoV eff N 2 {t) 




. (40) 


(30) 


V 2 (t)= f e l2X f^(t,x)dx, (31) 

Jo 

where 72 > 0 is a constant. 

Lemma 2. There exist positive constants 72 , ft, ft, 02 smc/z 
f/zflf f/ze following estimate holds for every solution to system 
© and © with © 


V 3 (t):= f e ' y3X f^ a; (t,x)dx, (41) 

Jo 

where 73 > 0 is a positive constant. 

Lemma 3. There exist positive constants 7,3 , S 3 . S 3 . 9 3 such 
that the following estimate holds for every solution to system 
© and © with © 

Vs(t) < - (ft + o(l))U 3 (f) - (ft + 0 ( 1 ))/^(ft 1 ) 

+ W(i) + / 2 (ft 1) + /£(ft 1)). (42) 

Proof of Lemma [3j Differentiating {41} gives, by {23}, {23} 
and {39} , that 


(32) 


V 3 (t) =BT 3 + / - 73 Q! p - 3a Pi e-^ x f 2 dx 


to L 


Proof of Lemma [2j Differentiating {4} and { 8 } with respect 
to t, gives that, 

(N(t) = k 1 -AP(t), Fi n [t) = k 2 ■ AP(t), 

[AP(f) = (i(t),JV(i))^{*) + §w\(l(t),N(t))^^- 


=BT 3 + (— 7 3 a Pe + o(l))U 3 (f), 


where 


PT 3 = - (e -73 a 

4p, 


+ ( tt Pe + °(l))/ PrT (C0). 


o(i))/L(u) 


(43) 


(44) 
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In order to estimate f Pxx (t, 0) or f Pxt (t, 0), essentially we 
need only to estimate Fi n (t ) and N(t), according to ( |33] and 
TO . On the other hand, © simply yields that 

N(t) = fciA P(i), F in {t ) = k 2 AP(t). (45) 



Therefore, from | 2 ], ( fTT] , ( |33] and ( [34] ). we have 

A P(t) = 0(l)(l(t) + N(t ) + f p (t, 1 ) + f pt (t, 1 )). (46) 

From ( |40] to ([46], we get 

f P , t (*, o) = o(i) (m+m+f p (t, i)+ f pt (t, i)). (47) 

Combining ([5]), | |20] , ( [35] ), ( |39] and ( |47] , we get further 

/ P „ (t, 0) =0(1) (/(f) + f p (t, 1) + f Px (f, 1)). (48) 

By ( [23] ) and ( [48] ). ( |44] becomes 

BT 3 =-{e-T*a Px +o{l))fi xm (t,l) 

+ 0 ( 1 ) (P(f) + />, 1 ) + fl(t, 1 )). (49) 

This implies the conclusion of Lemma [3] □ 

Step 4. The stabilization of l(t) and f p (t , •) in H 2 (0, 1). 
Finally, let the Lyapunov function be 

L(t) = V3 (f) + A3(V 2 (t) + A 2 (Vo(t) + AiVi(f))), (50) 

where A\ > 0 is such that © holds and A 2: Az > 0 
will be chosen later. Obviously, L[t) is equivalent to P(f) + 
\\U(t, 011 ^( 0 , !)■ Then - b y ©• ©• © and ( [50] , one can 
choose A 2 > 0 and A : > > 0 successively large enough such 
that 

L(t)^-(0 + o(l))L(t) (51) 

for some constant j3 > 0. We assume in a priori that 

for some small Eq > 0 such that |o(l)| ^ f in ©• Then 
L(t) ^ — §L(f), thus L(t ) ^ L( 0)e^^ t . Thanks to the 
assumption ( [13] , ( [52] can be satisfied for all t > 0 if e > 0 is 
small enough. The proof of Theorem [I] is thus complete. □ 

IV. Simulations 

Computing the time integration of the semi-discretized 
transport equations by finite volume with ODE45 routine of 
MATLAB, the stability result is achieved under the assump¬ 
tions of Theorems [l] Especially, the existence of k \, k 2 is 
guaranteed by a\ < 0 according to Proposition [2] The gain k\ 
is chosen to satisfy ( [IT] and ( [12] and k 2 is derived from the 
compatibility conditions ([9] and ( [TO] for which the inlet flow 
rate 73^(0) and the screw speed N( 0) with their respective 
time derivatives F in (Q) and N( 0) are computed with the help 
of the feedback ^ and the initial filling ratio f p ( 0 ). 

• Initial conditions: 

1° = 1.5 TO, f°(x) = 0.6905 + 0.025(1 - cos(irx)) + 
0.0117 sin( 7 ta;) 

• Setpoint values: 

l e = 1.37 TO, f pe = 0.6 




' 0 100 200 300 400 500 



915 °0 100 200 300 400 500 



Time (s) 

Fig. 4. Filling ratio || f p (t, ■) - / Pe || W 2 (01) 

• Gain values: 

h = 0.01, k 2 = 0.0001 

• Verification of assumptions ([TT]-(fT2] 

or = -0.0119 < 0 

V. Conclusion 

In this paper, we study the stabilization of a physical model 
for the extrusion process, which is described by conservation 
laws coupled through a dynamical interface. The exponen¬ 
tial stabilization is obtained for the closed-loop system with 
natural feedback controls through indirect measurements. The 
proof relies on Lyapunov approach. Numerical simulations 
are made as supplementary to the theoretical results. As 
a future work, it is interesting to study the controllability 
of boundary profile, i.e., to reach the desired moisture and 
temperature at the die under suitable controls. This problem 
is rather challenging for mathematical theory but also very 
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useful in applications. The proposed result is a first step 
towards controlling complex screw extrusion systems that 
include unknown parameters and time-varying disturbances 
acting on the pressure dynamics and the viscosity in the FFZ. 


Appendix 


Physical definition of the parameters 


L = 2 m 
B = 2.4 x 10” 6 m 4 
K d = 2.4 x 1CT 3 m 
C = 0.003 m 
g = 125 Pa s -1 
po = 350 kg m~ 3 
S eff = 0.014 m 2 
V eff = (S eff 


Extruder Length 
Geometric parameter 
Geometric parameter 
Screw Pitch 
Melt viscosity 
Melt density 
Effective area 
Effective volume 


Proof of Proposition [2] It suffices to consider the existence 
of ki satisfying Case 1) a i < 0. HD is true for ki — > 0 . 
Case 2) a i ^ 0 and 0261 0. If 62 = 0, HD is true either 
for fa —> Too or for fa —> —oo; If 62 7 ^ 0, HD is true 

either for fa -4 - \- or for ki —> -. Upon the existence 

o 2 o 2 

of fa, HI is true either for fa is arbitrary 1/0361 = 0 or for 
fa -T fpePoVefffa if 03&! 0. □ 
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